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We give the first self-consistent calculation of the effect of the scattered neutrino halo on flavor 
evolution in supernovae. Our example case is an 0-Ne-Mg core collapse supernova neutronization 
neutrino burst. We find that the addition of the halo neutrinos produces qualitative and quantitative 
changes in the final flavor states of neutrinos. We also find that the halo neutrinos produce a novel 
distortion of the neutrino flavor swap. Our results provide strong motivation for tackling the full 
multidimensional and composition-dependent aspects of this problem in the future. 
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I. INTRODUCTION 

In this paper we give the first self-consistent treatment 
of supernova neutrino flavor evolution that includes the 
"halo" [l] of neutrinos generated by direction-changing 
scattering in the supernova envelope. Neutrinos in the 
core collapse supernova environment are emitted primar- 
ily from the protoneutron star left by the collapse, which 
is considerably smaller than the envelope of the collapsing 
star. As neutrinos propagate outward, their flavor evolu- 
tion is determined, in part, by their coherent forward- 
scattering off of other neutrinos, almost all of which 
are emitted from the protoneutron star. However, the 
neutrino-neutrino coherent forward scattering contribu- 
tion to the potential that governs flavor transformation 
is sensitive to the intersection angle of interacting neu- 
trinos, such that the resultant potential is cx (1 — cos^), 
where is the angle between the incident neutrino trajec- 
tories. This proportionality has two direct consequence 
for the interactions of neutrinos in the supernova explo- 
sion. Neutrinos emerging from the protoneutron star 
will experience a significantly suppressed self-interaction 
once they have propagated more than a few protoneutron 
star radii. Neutrinos that have scattered at wide angles 
in the outer envelope of the explosion (halo neutrinos), 
while considerably less numerous than neutrinos emerg- 
ing form the core, experience no such suppression of the 
neutrino-neutrino forward scattering potential. This lack 
of geometric suppression leads to the result that during 
the first ~ 1 s of a core collapse supernova explosion, 
halo neutrinos can be the dominant source of neutrino- 
neutrino forward scattering potential in regions where 
active neutrino flavor transformation may take place [1]. 

Taken at face value, the halo changes the nature of 
flavor evolution calculations, converting them from ini- 
tial value problems into boundary value problems. The 
reason for this is that direction-changing scattering, in 
principle, can cause neutrino flavor information to prop- 
agate inward from a relatively large radius. Other studies 
of the halo have concentrated on stability of the neutrino 



flavor field in the accretion phase/shock reheting epoch 
of the supernova [2 , where composition, hydrodynamics- 
generated matter inhomogeneity, as well as inwardly 
propagating neutrinos currently preclude self-consistent 
calculations. In contrast, the very compact and cen- 
trally concentrated nature of the matter density distri- 
bution in 0-Ne-Mg core collapse supernovae makes the 
halo in these cases amenable to a self-consistent initial 
value treatment. Here we exploit this felicitous feature 
of 0-Ne-Mg core collapse and thereby take a step toward 
a more comprehensive treatment of neutrino transport. 

In fact, a longstanding and unresolved question in the 
physics of core collapse supernovae is the simultaneous 
and self-consistent solution of neutrino transport and fla- 
vor transformation. This subject has historically been 
approached by splitting the treatment of neutrinos in su- 
pernovae into two distinct limits: the Boltzmann trans- 
port limit, which contains all of the physical processes by 
which neutrinos are created, absorbed, and scattered in 
new directions; and the coherent forward-scattering limit, 
which governs the evolution of neutrino flavor states that 
are freely streaming. However, the existence of the neu- 
trino halo calls into question this separation that is at the 
heart of the current neutrino transport /flavor evolution 
paradigm. 

The halo shows that these two limits of neutrino trans- 
port confabulate in the general case. The neutrino halo 
effect itself is driven by the collusion of Boltzmann neu- 
trino transport together with the coherent forward scat- 
tering of neutrinos. The geometric structure of coher- 
ent neutrino-neutrino forward scattering strongly sup- 
presses the interaction of neutrinos that are propagat- 
ing on nearly co-linear trajectories, while enhancing the 
interaction of neutrinos that cross paths at large inter- 
section angles. Neutrino- nucleon/nucleus neutral cur- 
rent, isoenergetic interactions produce a population of 
neutrinos that may have scattered at wide angles, i.e. 
the halo neutrinos. This latter effect is dependent on 
the density and composition of matter in the envelope 
because of the coherent enhancement of neutrino scat- 
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tering on heavy nuclei [3^, W and, consequently, couples 
the halo neutrino population to the hydro dynamic and 
nuclear evolution of the supernova as a whole. Taken to- 
gether, coherent forward scattering and neutral current 
direction-changing scattering can give halo neutrinos dis- 
proportionate weight in determining the flavor evolution 
of all neutrinos ^Lj . 

At first glance, it may seem straightforward to include 
the halo neutrino population in calculations of neutrino 
flavor transformation, but the wide angle scattering that 
characterizes the trajectories of the halo neutrinos also 
makes their inclusion in such calculations difficult. Rel- 
ativistic Boltzmann transport of neutrinos has already 
been implemented in two and three spatial dimensions 
in simulations of core collapse supernovae [5]-[9]. This 
multi-dimensional transport treatment benefits from the 
fact that the relevant length scale for this kind of neutrino 
scattering is similar to the length scale relevant to resolv- 
ing the hydrodynamic evolution of the envelope. Unlike 
Boltzmann neutrino transport, the neutrino coherent for- 
ward scattering limit demands that the complex phases 
of neutrino flavor state high frequency oscillations be fol- 
lowed on length scales as small as ~ (a few) cm. The 
very short length scales associated with self-consistently 
following neutrino flavor evolution have stymied efforts to 
expand simulations to multiple spatial dimensions. Cur- 
rent state of the art calculations in this field are pred- 
icated on the assumption that the only neutrino states 
that need to be followed are outwardly-directed. The dis- 
covery of the neutrino halo has shown this assumption 
to be untrustworthy during the core collapse explosion 
epoch, and appropriate only during the post-explosion 
neutrino driven wind epoch. 

This does not mean, however, that the effects of the 
neutrino halo on flavor transformation are outside the 
realm of consideration. The low mass of the 0-Ne-Mg 
core collapse supernova progenitors results in a prompt 
and spherically symmetric explosion. This allows us to 
use our extant, spherically symmetric and multi- angle 
flavor transformation calculations for this particular case 
at early times. Further, a calculation of neutrino flavor 
evolution can be made so long as there are few enough 
halo neutrinos on inwardly-directed trajectories that they 
do not contribute significantly to forward scattering po- 
tentials (as is the case during the neutrino driven wind 
epoch) [T. Finally, there must be some reasonable ex- 
pectation that all neutrino flavor states can be specified 
from the start of the calculation. 

The neutronization burst stage of the 0-Ne-Mg core 
collapse supernova satisfies these criteria because of its 
centrally condensed matter envelope. Figure [l] shows the 
density versus radius profile for the 0-Ne-Mg configura- 
tion of Reference [10^ 11 . The region of the envelope that 
provides significant numbers of neutrinos to the halo pop- 
ulation is located within a radius of ~ 1000 km [1]. At 
the edge of this region the matter density drops precip- 
itously, and ceases to scatter enough neutrinos into the 
halo to make a significant contribution to the forward 
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FIG. 1: The matter density profile taken from Reference [10] . 
The profile is quite centrally concentrated, so much so that the 
scattered halo is only populated with an appreciable number 
of neutrinos inside a radius of ^ 1000 km. 



scattering potentials. A calculation of flavor transfor- 
mation that begins outside this point may safely neglect 
halo neutrinos on inwardly-directed trajectories. 

Previous studies have shown that flavor transformation 
in the 0-Ne-Mg neutronization burst does not proceed 
until a radius of ~ 1100km p!2UT4) . This is due to the 
large flux of neutrinos emitted during the neutronization 
burst. What this means for this particular case, where 
no flavor transformation has taken place inside the re- 
gion which scatters neutrinos into the halo, is that the 
flavor states of all of the halo neutrinos can be deter- 
mined from the flavor states of neutrinos emerging from 
the core. This allows us to perform our spherically sym- 
metric, multi- angle, flavor transformation calculations by 
starting at a radius where backwards going neutrinos are 
negligible. 

In what follows we discuss the methodology of our cal- 
culation of in Section II, and results in Section III. We 
give conclusions in Section IV. 



II. METHODOLOGY 

For the neutrino emission during the neutronization 
neutrino burst, we use the results of Reference [10 to 
set the spectral energy distribution of all flavors of neu- 
trinos. The neutrino emission parameters for each time 
slice we consider are show in table |lj We use these pa- 
rameters to fit the total neutrino emission to a normal- 
ized Fermi-Dirac spectrum for each flavor of neutrino and 
anti-neutrino 
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TABLE I: Neutrino emission parameters for the initial spectra 
used in our calculations. 

where we take rjj^ = 3 and Tj^ = F2{r]j^) (Ejy) / Fs{r]j^) . Here 

poo 

FJri) = / z ^ dx. (2) 

"^^^ Jo exp(x-7?) + l ^ ^ 

For the purposes of this study we have chosen the follow- 
ing neutrino mixing parameters: neutrino mass squared 
differences Am| = 7.6 x 10"^ eV^ and Aml^^ = 2.4 x 
10"^ eV^; vacuum mixing angles O12 = 0.59, O23 = 7r/4, 
6is = 0.152; and CP- violating phase S = 0. 

Fundamentally, the calculation that we are perform- 
ing in this case is an initial value problem specified at 
a fixed initial radius. To begin, we specify a unique in- 
tensity and spectral energy distribution for each flavor of 
neutrino and anti-neutrino on a surface of fixed radius, 
and simultaneously solve the non-linearly coupled equa- 
tions of motion for the evolution of neutrino flavor states 
in the coherent forward scattering limit. The require- 
ment of simultaneity in the calculation of neutrino flavor 
states forces the entire solution to move outward, in lock- 
step along the radial coordinate, for all neutrinos. By 
way of contrast, the neutrino halo is a manifestly multi- 
dimensional phenomenon, with neutrinos moving in all 
directions. Including the halo in an initial value problem 
formalism such as ours raises a number of thorny issues 
that must be addressed. 

The most stringent requirement is that the calculation 
must be limited entirely to the region where halo neu- 
trinos propagating inward, opposite to the direction of 
the calculation itself, are truly negligible. So that we 
may compare the effect of the halo neutrinos to previ- 
ous work, we employ the density and composition pro- 
file found in Reference [10], which is fit to the collapse 
progenitor profiles of Refs. [151 IS] when flavor transfor- 
mation calculations extend past the outer radial coor- 
dinate of the profile of Reference [10 . To compute the 
neutral current neutrino scattering off of nucleons and 
nuclei we again use the energy dependent cross sections 
found by Reference [4 . For any scattering that takes 
place outside the simulation volume of Reference [lOj, 
and hence originates in the density profile of the progen- 
itor p!5^, , we assume that the composition of the enve- 
lope is entirely helium until the hydrogen burning shell 
is reached at r = 1090 km, and pure hydrogen outside 
of this point. In this calculation we neglect direction- 
changing neutrino-electron scattering. While the cross 
sections for neutirno-electron scattering are comparable 
to the neutrino-neutral current scattering cross sections 
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FIG. 2: The neutrinosphere radius, as defined by the surface 
of optical depth r = 1, shown as a function of neutrino energy. 



for helium and free nucleons, the low mass of electrons 
relative to the energy of neutrinos in the neutronization 
neutrino burst produces scattering that is strongly for- 
ward peaked. For this reason, neutrinos that have under- 
gone electron scattering do not contribute significantly to 
either the halo neutrino population or potential. 

In keeping with our previously stated 1 % criterion P, 
we find that outside a radius of 850 km the neutrinos 
scattered onto inward trajectories do not contribute more 
than 1 % to the magnitude of the neutrino self-coupling 
potential. We use this result to define a "halosphere" , 
at the radius Rn = 850 km for this model, which is the 
surface outside of which the propagation of halo neutrinos 
may be taken to be in the outward direction without 
impacting the dynamics of flavor transformation. 

A second requirement that must be met is that we must 
make a physically motivated choice for the initial flavor 
states of the neutrinos that have been scattered into the 
halo population. As mentioned in Section I, multi-angle 
suppression prevents the onset of neutrino flavor trans- 
formation. We find that the neutrino driven multi-angle 
suppression alone is sufficient to suppress collective oscil- 
lation out to a radius of r = 1100 km, and that this figure 
is little changed by the addition of the halo to our cal- 
culations. Therefore, no neutrinos within this radius will 
have had the opportunity to engage in collective flavor 
oscillation. Because this radius is outside the halosphere 
surface, all neutrinos scattering into the halo within the 
halosphere will be in their original flavor state. Further- 
more, because the halo is populated by neutral current 
processes that are flavor blind, the halo neutrinos will 
not change flavor after scattering. This uniquely deter- 
mines the flavor states of all neutrinos, whether emitted 
directly from the core or scattered in the halo, at the 
surface of the halosphere. 

In our previous papers on this subject we have used 
the convention that all neutrinos are emitted isotropi- 
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FIG. 3: Initially neutrinos are emitted isotropically from 
energy- determined neutrinosphere surfaces, R^"^ . After emis- 
sion, some neutrinos are scattered at wide angles into the 
neutrino halo (e.g. v\^i jvv scattered to location at radius r 
and angle of incidence ^k), while the remainder continue on 
unimpeded to the halosphere surface at radius i?H. The tra- 
jectories of all neutrinos emerging from the halosphere surface 
are characterized by the angle they make relative to the 
outward unit normal on this surface and their propagation 
distance /k from the halosphere surface. 



cally from a hard spherical shell called the neutrinosphere 
(which is just above the surface of the proto-neutron 
star). A primary criticism of this picture has been that 
the radius of the neutrinosphere surface itself is depen- 
dent on neutrino energy, particularly for high energy neu- 
trinos, and the emission region is extended for these neu- 
trinos. Adapting our calculation to accommodate dif- 
fering neutrino emission spectra along different emission 
trajectories has afforded us the opportunity to rectify this 
shortcoming. 

A neutrinosphere surface can be defined where the op- 
tical depth against scattering is r = 1. For a typical neu- 
trino (average energy) with our density profile this cri- 
terion gives a radius = 60 km. Each neutrino energy 
can be assigned an appropriate "neutrinosphere" where 
the optical depth is unity. Shown in Figure |2] is the com- 
puted radius of each r = 1 neutrinosphere surface for 
neutrinos of a given energy. 

For each neutrino energy we compute the neutral cur- 
rent scattering into the neutrino halo at the surface of 
the appropriate, energy dependent, neutrinosphere. This 
means that the halo region itself is not sharply defined, as 
low energy neutrinos are scattering into the halo at radii 
where high energy neutrinos are still inside of their re- 
spective neutrinospheres. The intensity of neutrino emis- 
sion on the surface of the halosphere is specified by the 
neutrinos emerging through the surface along a given tra- 
jectory with angle i^k relative to the outward unit normal. 
The ensemble of neutrinos at the halosphere surface is 
populated by emission from the neutrinospheres as well 
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FIG. 4: The neutrino emission intensity as a function of im- 
pact parameter at the halosphere surface. Line style indicates 
different neutrino energies: lOMeV (solid), 30MeV (dash- 
dotted), 50MeV (dashed), 70MeV (dotted). 



as from the neutrino halo. Figure [3] shows a cartoon rep- 
resentation of the physical set up of the initial conditions 
where i^k is the neutrino emission angle relative to the 
surface of the halosphere, r is the distance from the cen- 
ter of the supernova, is the angle of intersection of the 
neutrino trajectory with the outward unit normal, and /k 
is the propagation distance of vy^ from the surface of the 
halosphere to the location at radius r. 

This approach carries the explicit assumption that 
there is a clear break between the multiple scattering 
regime, where the optical depth of neutrinos is r > 1, 
and the single scattering regime, where r < 1, at the 
(energy dependent) neutrinosphere surface. We take the 
emission of neutrinos of a given energy to be isotropic 
at the neutrinosphere surface (which is the limit in the 
case of multiple neutrino scattering), and employ only 
single scattering to compute the angular dependence of 
the halo neutrino population outside the appropriate neu- 
trinosphere. While this approach is somewhat crude, it 
produces results that are in remarkably good agreement 
with sophisticated models of neutrino transport in super- 
novae gllSHZl. 

The intensity of the neutrino emission at the halo- 
sphere surface is shown in Figure |4] Note that because of 
the existence of multiple neutrinospheres and the diffuse 
scattering of halo neutrinos, we have chosen to param- 
eterize the individual neutrino trajectories by their im- 
pact parameter, b = Rusini^]^^ relative to the center of 
the proto-neutron star. For a given neutrino energy Ejy, 
the neutrino emission "intensity", 1^ {Ejj)^ shown in Fig- 
ure |4] is defined in this case to be related to the neutrino 
number density, {Ejy)^ outside the halo sphere surface 
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by the relation, 



"2 dcosO]^. (3) 



III. RESULTS 



Our results show that there are both quantitative and 
qualitative changes in the structure of the neutrino fla- 
vor transformation in the presence neutrino halo. We 
consider two different time slices of the neutronization 
burst, one 7 ms post core bounce at the height of the 
Uq luminosity, and another 15 ms post bounce when the 
fluxes of other species of neutrino have begun to rise ap- 
preciably. All neutrino emission parameters are taken 
directly from Reference [10 . 

The salient question raised by the addition of halo 
scattering to our flavor transformation calculations is 
whether or not this effect has any observable conse- 
quence. To address this, we have taken the fluxes of 
neutrinos generated by our flavor transformation calcu- 
lations for both the halo and no halo scattering cases and 
used the SNOwGLoBES software package [17 to model 
the detected signal corresponding to the two spectra seen 
in Figure [5| The primary signal for the neutronization 
burst epoch will be in the v sector, so we have chosen 
to use a liquid Argon model detector that is most sensi- 
tive to the Vq flux via charged current capture [18 . The 
emission angle averaged results of our calculations for the 
15 ms post bounce time slice are shown in Figure [5] The 
results in Figure [5] show that the addition of the halo 
scattering has produced a swapped population of mass 
state 2 neutrinos above an energy of 20MeV. Because 
mass state 2 has a larger electron flavor component than 
mass state 3, the swapping of neutrinos into mass state 2 
shown if Figure [5] is expected to produce a signal which 
should be observable. Shown in Figure [6] is the difference 
between the modeled detector signals with and without 
halo neutrino scattering, for the portion of the neutron- 
ization burst where the emitted fluxes are similar to what 
was used in the flavor transformation calculations shown 
in Figure |5j a window of ~ 20 ms. Figure [6] demonstrates 
that the effect of the halo scattering has been to produce 
a clear swap feature that is detectable in the received 
supernova neutronization burst signal. 

Figures 7-10 show in detail the results of our flavor 
transformation calculations. These figure show the prob- 
abilities for a neutrino or anti- neutrino, initially in the 
electron flavor state, to occupy each of the neutrino mass 
basis states at large radius, after neutrino flavor trans- 
formation is complete. We define this probability to be: 



K^^i (initial) I (^final))!^ 



(4) 



Here l^^y. (rinitiai)) is the wave function for a neutrino of 
flavor i = e, /i, r at the radius rinitiai where the neutrino 
flavor states are initialized, and 1^^^^ (rfinai)) is the wave 



function for a neutrino mass state a = 1, 2, 3 at the ra- 
dius rfinai where the calculation ends. Figure [7| shows the 
mass state occupation probabilities for electron neutrinos 
during the 7 ms time slice, comparing the cases where 
the initial states of all neutrinos are prepared in with 
and without scattering into the halo population. Fig- 
ure [8] shows the mass state occupation probabilities for 
electron anti- neutrinos during the 7 ms time slice, com- 
paring the cases where the initial states of all neutrinos 
are prepared in with and without scattering into the halo 
population. Similarly, Figures |9] and [lO] display the com- 
parison of the Halo and No Halo cases for the i^e/^e mass 
state occupation probabilities, respectively, for the 15 ms 
time snapshot. 

For the cases where no halo scattering is included, 
the emission from the surface of the neutrinosphere is 
isotropic, and the trajectory bins that lie in the halo re- 
gion are unpopulated. For each time snapshot the Halo 
vs. No Halo calculations employ identical energy spectra 
for all neutrino flavors, and identical binning schemes for 
both emission angle and neutrino energy. Furthermore, 
the emission trajectory binning for both calculations was 
chosen so that there would be an exact match to the bins 
on the surface of the neutrinosphere employed in previ- 
ous calculations [19] , with an equal number of additional 
trajectory bins added to accommodate the halo neutri- 
nos. 

The most noticeable feature of Figures [7] and [8] is that 
the mass state 3/2 swap is emission angle independent 
when neglecting the halo effect, and very much depen- 
dent on the neutrino trajectory when the halo is included. 
This is a surprising result, as the analytic work which has 
studied the formation of these flavor swaps does not find 
any dependence of the swap energy on the emission an- 
gle. Fascinatingly, with the halo the swap interface is 
pushed through Ey = (X) dX large emission angles, and 
reappears in the anti- neutrino sector, creating a swap in 
the anti-neutrinos which is not present in the absence of 
the halo. Even more interesting is that unlike the mass 
state 3/2 swap, the swap between mass state 2/1 exhibits 
an apparent lack of emission trajectory dependence. 

For the Figures [9] and 10 the disparity of between num- 
ber fluxes of Vq and Vq has diminished after the burst 
earlier reached peak luminosity,which is typical of the 
neutronization burst of core collapse models from which 
these spectra are drawn As can be seen in Fig- 

ure [To) this shift in the spectral properties of the emis- 
sion has moved the bulk of the spectral swaps into the v 
sector when scattering into the halo is neglected. How- 
ever, when the effect of the halo is included the spectral 
distortion it produces shifts the mass state 3/2 swap back 
into the v sector. This creates a population of swapped 
mass state 2 neutrinos above ~ 20MeV, which was not 
present when the calculation was performed in the ab- 
sence of the halo. As we discussed earlier in this section, 
this distortion of the spectral swap created by the halo 
scattering produces a clearly detectable swap signal in 
the V sector for this time snapshot. 
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FIG. 5: A comparison of the emission angle averaged results of flavor transformation calculations with the halo neutrinos 
included and with halo scattering neglected. Left panel: the calculation including the halo, mass basis (key top right, inset) 
neutrino energy distribution functions versus neutrino energy. The dashed curve gives the initial v energy spectrum. Right 
panel: the calculation neglecting halo scattering, mass basis (key top right, inset) neutrino energy distribution functions versus 
neutrino energy. The dashed curve gives the initial v energy spectrum. Both panels show the final state of neutrino fiavor 
transformation at a radius of r = 12000 km. 
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FIG. 6: A comparison of the modeled event rate for detected 
captures in a 17kt liquid Argon detector between calcu- 
lations with and without the scattered neutrino halo. The 
spectral distortions created by the halo produce a clear swap 
signature between 20 — 30 MeV, which constitute ^ 15 ad- 
ditional Ve events in this 20 ms time slice of the supernova 
signal. 



IV. THEORY 

The spectral distortions found in our calculations raise 
a question: Do the halo neutrinos, though few in number, 
nevertheless alter the qualitative and quantitative char- 
acter of collective neutrino oscillations? The answer: At 
7 ms in our model the halo primarily affects the collec- 
tive oscillations of neutrinos propagating at large impact 
parameters; but 8 ms later the halo neutrinos completely 



re-determine the course of neutrino flavor oscillation for 
all emission trajectories. This result underscores the ne- 
cessity for a self-consistent numerical treatment of this 
nonlinear system. 

The twisting of one of the swap surfaces through the 
trajectory space has several direct consequences. The 
first is the shift in the swap energies. When the halo 
effect is included in the 7 ms post bounce case, a high en- 
ergy tail of remains unswapped in the neutrino sector. 



Figure [TT] shows this feature in the total angle- aver aged 
energy spectra for electron neutrinos projected into the 
three mass states for our simulation with and without the 
halo. The total number of neutrinos in each mass state 
for both the halo and no-halo cases are nearly identical 
(there are small differences on the order of ~ 0.1 %, owing 
to slight increases in the adiabaticity of flavor evolution 
when the halo is included). With the halo the number 
of neutrinos that remain in mass state 3 at high energy 
causes the swap between i'2,/'^2 to form at lower energy. 
Consequently, this also lowers the swap energy for mass 
states 

This can be understood simply from the equations of 
motion. The collective flavor oscillation which creates the 
swaps (called the Regular Precession Mode) in this ex- 
ample posseses two conserved constants of the motion, 
effective lepton numbers for each of the mass-squared 
splittings [20 . Because the scattering of neutrinos into 
the halo does not change the spectral shape of the entire 
ensemble of neutrinos, one might reasonably expect that 
the conserved lepton numbers that describe the flavor 
evolution of the neutrinos to remain unchanged by the 
presence of the halo. Indeed, this is what is found in our 
calculations. Following the convention of Reference [20] . 
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FIG. 7: Results of flavor transformation calculations with the halo neutrinos included. Initial emission spectra are taken 
from [To], 7ms post core bounce. Left panels: electron neutrino occupation probability Pu^u^ (color/shading key at top of 
panel), where x = 1, 2, 3 is the neutrino mass eigenstate, shown as a function of impact parameter, b in units of neutrinosphere 
radius Ri^ = 60 km, and neutrino energy, E in MeV, plotted at a radius of r = 12000 km. Right panels: electron anti-neutrino 
occupation probability Pu^u^ (color/shading key at top of panel), where x = 1, 2, 3 is the anti-neutrino mass eigenstate, shown 
as a function of impact parameter, b = Ru sini^k, and neutrino energy, E in MeV, plotted at a radius of r = 12000 km. 
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FIG. 8: Same as Figure [7| for z^e's. 



the conserved lepton numbers and L3 for the atmo- 
spheric and solar mass squared splittings, respectively, 
are identical for calculations with and without the pres- 
ence of the halo. As a consequence of this conservation 
law, the presence of additional neutrinos in mass state 



three at high energy necessitates a concomitant reduc- 
tion in the swap energy seen in Figure [7[ 

The conserved lepton numbers are also nearly identical 
for the calculations performed for the neutrino emission 
15 ms post core bounce. Although the total occupation of 
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FIG. 9: Same as Figure|9]for z^e's during the 15 ms time snapshot. 



mass state 2 in the neutrino sector is manifestly different in the v sector can significantly alter the v contributions 

for the Halo vs. the No Halo case, the contributions of to and L3. To satisfy the conservation of lepton num- 

z/'s and z>'s to the magnitude of and L3 have a relative ber demanded by the equations of motion, the spectral 

sign difference. Once the number fluxes of v begin to rise distortions in the v sector must shift an appreciable num- 

during the neutronization burst, the spectral distortions ber neutrinos into the appropriate mass Eigen states to 
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FIG. 10: Same as Figure[9]for z^e's during the 15ms time snapshot. 



balance the spectral distortions in the v sector. 

The agreement on the conserved lepton numbers be- 
tween the Halo and No Halo calculations strongly sug- 
gests that the distortion of the swap surface seen in Fig- 
ures 7 - 10 is due to differences in the geometry and 



spatial distribution of neutrinos between the halo and 
no halo calculations, and is not a product of neutrinos 
following a different equation of motion through flavor 
space in the presence of the halo. The Regular Pre- 
cession Mode is characterized by strong non-linear cou- 
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FIG. 11: Emission angle averaged results of flavor transformation calculations at 7 ms post core bounce comparing calculations 
with and without halo neutrino scattering. Left panel: mass basis (key top right, inset) neutrino energy distribution functions 
versus neutrino energy when halo scattering is included in the calculation. The dashed curve gives the initial Ue energy 
spectrum. Right panel: mass basis (key top right, inset) neutrino energy distribution functions versus neutrino energy when 
halo scattering is excluded from the calculation. The dashed curve gives the initial i>e energy spectrum. Both panels show the 
final state of neutrino flavor transformation at a radius of r = 12000 km. 



pling of neutrino flavor states, leading to all neutrinos in 
the ensemble oscillating in flavor space with the same 
frequency, cjpr. In the single angle approximation all 
neutrinos either align or anti-align with a mass eigen- 
state as the neutrino self-coupling potential decreases, 
and cjpr = cjswap = Am^/2^swap- The sense of the swap 
depends on the original alignment of the neutrino mass 
states and whether the individual vacuum oscillation fre- 
quencies of those neutrinos are greater or less than cjgwap- 
In the calculation presented here, which employs the nor- 
mal neutrino mass hierarchy, a neutrino with vacuum os- 
cillation frequency uy > Uswap will remain in the heavy, 
z/3, mass eigenstate. 

Requiring that all neutrinos in the ensemble, across 
all emission trajectories, agree on cjpr, as suggested by 
our results when the halo is included, admits an interest- 
ing solution. In the single angle formulation, a neutrino 
which is precisely at the swap energy satisfies the crite- 
rion cjpr — oJy = 0. In the region of the envelope where 
the swap is forming, e.g., where {\Hi^j^\) ~ cjpr, there 
may be a significant difference between a neutrino's vac- 
uum oscillation frequency and its instantaneous flavor 
oscillation frequency. The latter arrises from dispersion 
in the forward scattering potentials along different emis- 
sion trajectories. This suggests a natural modification of 
the swap criterion which accommodates the additional 
dispersion present when the multi-dimensionality of the 
supernova environment is accounted for: 



dent swap frequency. 



(CJV + ^H^u + Ai^e + AHy) = 



(5) 



This is directly equivalent to a shift in the swap frequency 
for each emission trajectory. By grouping the dispersion 
terms with cjpr, we define an effective, trajectory depen- 



{AH,, + AH^ + AHy) 



(6) 



This explains neatly the effect seen in Figures 7 - 
10. In Figure [7| as the dispersion effect increases with 
increasing impact parameter, the mass state 3/2 swap 
energy is pushed to much higher values, indeed cross- 
ing into the anti-neutrino sector, which corresponds to 



^swap < 0. Figure 12 shows the isocontour which satis- 
fies cDswap — = Oior the mass state 3/2 swap. This 
is plotted beneath the mass state occupation probability 
data. 

When the isocontours for both swaps are considered, 
we recover the shape of the mass state 2 occupation prob- 
ability which has produced the detectable swap feature 
for the 15 ms snapshot, shown in Figures [9]and[l0l In 
Figure 13 both the isocontrous for the mass state 3/2 
and 2/1 swaps are shown. This is plotted alongside the 
mass state occupation probability data for mass state 2. 
Because the mass state 2/1 swap is formed by neutrino 
oscillations in the AttIq sector, the swap itself is formed 
when the neutrino self-coupling potential has a smaller 
magnitude than when the mass state 3/2 swap is formed. 
Consequently, the mass state 2/1 swap is formed much 
further out in the envelope of the supernova. At these 
larger radii, the dispersion effects which distort the swap 
surfaces are reduced, resulting a flavor swap that is less 
sensitive to the neutrino trajectory. 
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FIG. 12: Top panels: the structure of the mass state 3/2 swap in the neutrino (left) and ant i- neutrino (right) sectors, displayed 
in terms of the mass state 3 occupation probability for the 7 ms post bounce calculation (notation and axes as in Figure y\ . 
Bottom panels: shown in black is the contour which satisfies the emission trajectory dependent swap criterion of Equation [5| 
for neutrinos (left) and anti- neutrinos (right). The thicker, lighter color contour shows the locations where a or has a 
50% probability to occupy mass state 3. The contours displayed are selected at a radius where \Hvv\ ^ cjpr is satisfied for a 
given impact parameter, h. 



V. CONCLUSION 

We have made the first mult i- angle calculation of 
neutrino flavor evolution in the supernova environment 
which includes the population of neutrinos scattered into 
the diffuse neutrino halo. We have shown that there 
are qualitative differences in the neutrino oscillation pat- 
terns as a function of angle and energy. These qualita- 
tive changes may also have potentially detectable conse- 
quences for a received neutronization burst signal. 

This calculation was made possible by a confluence 
of physical circumstances present during the neutroniza- 
tion burst of an 0-Ne-Mg core collapse supernovae. The 
combination of multi-angle suppression of neutrino flavor 
transformation deep within the envelope, and the precip- 
itous drop in matter density just outside of the volume 



where multi-angle suppression ceases to operate, creates 
a unique situation where the halo neutrinos moving on 
inward directed trajectories are negligible at the radii 
where neutrino flavor transformation takes place. This 
configuration allows for the inclusion of the neutrino halo 
within the present, initial value problem framework of 
computational models of supernova neutrino flavor trans- 
formation. 

The results of our calculations show that the neutrino 
flavor swap, the clearest signature of neutrino collective 
oscillations, is a phenomenon dependent on the geometry 
of the neutrino flavor transformation environment, and 
hence of the geometry of the supernova envelope. The 
geometric dependence of the swap energy is a feature of 
neutrino collective oscillations which has been completely 
overlooked by previous calculations. The dependence of 
the halo neutrino fiux on both the matter density and 
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FIG. 13: Top panels: the structure of the mass state 2 occupation probabihty in the neutrino (left) and anti- neutrino (right) 
sectors for the 15 ms post bounce calculation (notation and axes as in Figure [?]). Bottom panels: shown in solid black is the 
contour which satisfies the emission trajectory dependent swap criterion of Equation [s] for the mass state 3/2 swap for neutrinos 
(left) and anti-neutrinos (right). Shown in dashed black is the contour which satisfies the emission trajectory dependent swap 
criterion of Equation [5] for the mass state 2/1 swap for anti-neutrinos (right). The thicker, lighter color contour shows the 
locations where a Ve/ve has a 50% probability to occupy mass state 2. The mass state 3/2 contour and the mass state 2/1 
contour are selected at differing radius so that for both contours the condition \Huu\ ^ cjpr is satisfied individually for each 
combination of precession frequency, cjpr, and impact parameter, b. 



composition of the envelope now implies that the shape 
and spectral properties of the flavor swaps created in the 
explosion may also bear a dependence on the envelope as 
well. 

Out results demonstrate the necessity of a self- 
consistent numerical approach in modeling collective os- 
cillations in the fiercely nonlinear environment of stellar 
collapse. For example, we found that while the small 
number of halo neutrinos has little effect on collective 
oscillations at 7 ms, only 8 ms later these halo neutrinos 
significantly alter both the qualitative and quantitative 
process of flavor oscillation and swap formation. 

The spectral distortions that the halo creates have 
potentially detectable consequences for a received neu- 
tronization neutrino burst signal here on Earth. The 



halo shifts the apparent swap energies of the neutrino 
signal, sometimes across the boundary between v and 
v sectors. This is accomplished without changing the 
overall spectral properties of the initial neutrino states. 
Attempts |T9l [21] to reverse engineer a swap signal to 
extract information on the supernova environment must 
take account of the halo effect. 

Finally this work gives a tantalizing glimpse of the 
new phenomenology which has emerged in spherical sym- 
metry. The Halo is, of course, a fundamentally multi- 
dimensional phenomena and these results strongly mo- 
tivate attempts to expand the dimensionality of neu- 
trino flavor transformation calculations. Further, there 
are a multitude of different progenitor models which pro- 
duce distinct signals [22 during the supernova explosion. 



14 



While we have presented an example of a single case here, 
a general solution for the effect of the scattered halo 
on flavor transformation in the explosion is intractable 
at present for later times and more massive supernova 
progenitors. The distortion of the swap energy surface 
through emission angle space is a phenomenon that re- 
veals how robust the collective oscillation modes of neu- 
trinos can be. Further, it exposes as false the fundamen- 
tal assumption of the single-angle approximation: that 
individual neutrinos of the same energy and initial fla- 
vor state, following the same equation of motion, should 
be in identical flavor states. If the distortion of spec- 
tral swaps persists in the presence of large numbers of 
inward directed neutrinos, it will be exacerbated by local 
sources of dispersion in the neutrino self-coupling poten- 
tial, such as halo neutrino reflections off of turbulence 
driven cold matter accretion plumes. This work raises 
the point that the coherent forward scattering of neu- 
trinos and the Boltzmann transport of neutrinos do not 
belong in the separate camps to which they have appor- 



tioned in the supernova environment. The inclusion of 
the neutrino halo in flavor transformation calculations is 
a zeroth order attempt to bridge this gap, and the re- 
sults point directly to the need for full quantum kinetic 
treatment for the general case of neutrino transport in 
super novae. 
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